Collision of One-Dimensional Nonlinear Chains 
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We investigate one-dimensional collisions of unharmonic chains and a rigid wall. We find that the 
coefficient of restitution (COR) is strongly dependent on the velocity of colliding chains and has a 
minimum value at a certain velocity. The relationship between COR and collision velocity is derived 
for low-velocity collisions using perturbation methods. We found that the velocity dependence is 
characterized by the exponent of the lowest unharmonic term of interparticle potential energy. 
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I. INTRODUCTION 

In collisions between two bodies which have internal 
degrees of freedom, some of the initial translational en- 
ergy is transformed into internal energy of the two bod- 
ies. This is the major cause of energy dissipation. To 
characterize the macroscopic features, a phenomenologi- 
cal parameter, the coefficient of restitution (COR), 



K r 
~K~i 



(1) 



where Ki and K r are translational kinetic energy before 
and after the collision, respectively, is commonly used. 
Recent studies of collisions are mainly focused on the 
determination of r\ from microscopic mechanisms. 

Hertz developed the theory of collision between fric- 
tionless elastic bodies jlj || based on his static theory of 
elastic contact ||. In the theory, it is assumed that in 
low- velocity collisions, the deformation of colliding bodies 
is given by the static theory and the production of vibra- 
tion is totally ignored. Hence, the theory gives no infor- 
mation on energy dissipation. Plastic deformation is one 
of the possible means of kinetic energy dissipation during 
collisions. Taking this into account, the energy dissipa- 
tion rate 1 — r\ is found to increase with collision velocity 
with a power law with the exponent 1/2 [|J. Considering 
viscoelastic properties, the dissipation rate increases with 
the exponent 2/5 5, 0. These results were confirmed 
experimentally ||, £, Rj. 

The results presented above are based on quasi-static 
approximation, hence it is expected that they are re- 
stricted to low-velocity collisions. We believe that more 
general results will be obtained through microscopic sim- 
ulations (nj |n| [[3|. 

Sugiyama and Sasaki |Q and Basilc and Dumont 
p5[ considered collisions between simple one-dimensional 
chains and a rigid wall. A chain is composed of n identical 
point particles which interact with nearest-neighbor par- 
ticles. If we choose linear force as the interaction between 
the particles, the COR is independent of collision velocity 
and approaches unity in the thermodynamic limit. 

For collision between metallic bodies, plastic deforma- 
tion usually occurs even if collision velocity is low. The 
deformation exceeds the elastic regime in which Hooke's 



law remains valid |J, [lfj . Hence nonlinear effects are im- 
portant in realistic collisions. In this study, we choose 
nonlinear force as the interparticle interaction and per- 
form the one-dimensional simulations of collision between 
a "nonlinear chain" and a rigid wall. We find that the 
COR of this collision has a minimum value at a certain 
velocity and derive the velocity dependence of the COR 
for low- velocity collisions using perturbation methods. 

This paper is organized as follows. In the next section, 
we briefly review the studies of Sugiyama and Sasaki 
and Basile and Dumont Jl5| . In section 3, we present our 
results for collision between nonlinear chains and a wall. 
Finally, we summarize our results. 



II. COLLISIONS OF ONE-DIMENSIONAL 
HARMONIC CHAINS 

First we briefly discuss the collision of one-dimensional 
harmonic chains (model A) with a rigid wall as discussed 
by Sugiyama and Sasaki |14j . Consider a chain composed 
of n identical point particles labeled j = 1, 2, • • • n. Each 
particle in the chain is linked to nearest-neighbor parti- 
cles with a Hooke's spring, as illustrated in Fig.[l]. The 
Hamiltonian of this system is written as 



3 = 1 3 = 1 



(2) 
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FIG. 1: Schematic of a one-dimensional chain and a rigid wall. 
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where Xj and Xj are the position and velocity of the j- 
th particle, respectively, and £ is the natural length of 
the springs. We assume the chain is homogeneous and 
all springs have the same spring constant k(= maj 2 ). V w 
represents the hard-core potential of a rigid wall located 
at x — where collision takes place. During the collision 
process, only the particle j = 1, which is at the end of the 
chain, interacts with the wall since the positional order 
of the particles is always kept as x n < x„_i < • • • < X2 < 
X\. In the following discussion, we assume that the wall 
is so rigid that the particle X\ simply reverses its velocity 
as — x\ — ► x\. 

Although there is no dissipation term in the Hamilto- 
nian, vibration energy which remains after the collision 
is regarded as the "dissipated" portion of energy. The 
COR 7/ is therefore evaluated as 



V 



1 



E, 



vib. 

~E~ 



(3) 



where E is the total energy and E v n>. is the vibration 
energy which remains after the collision. 

In the absence of V w , the equation of motion is written 

as 



x = -w 2 (Ax + b), 

where x = (x\, X2, ■ ■ ■ , x n ), b = {—I, 0, 
is an n x n matrix of the form 



A = 



/ 1 -1 
12-10 
12-10 
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(4) 
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(5) 
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Taking a principal-axis coordinate system in which A is 
diagonal, internal vibration of the chain can be repre- 
sented using n nonintcracting fundamental modes. In 
the presence of V w , collision between the chain and the 
wall is realized as follows. Assume that the particle 
j = 1 collides with the wall f(n) times at t%,t2, ■ ■ ■ , tf(n)- 
The fundamental modes describe equi-energy elliptical 
orbits in phase space. The orbits are discontinuous at 



t = t%, ti, • 



■ t 



/(")• 



For the following numerical simulations, we choose ini- 
tial conditions 



-wo (j 



1,V 



(6) 



Before collision, the chain has no internal vibration, i.e., 
zero temperature. Figure || shows the collision between 
the wall and the chain with n = 19. Each line is a trace of 
the trajectory of a particle. In the plot, units on the time 
axis are taken as r = (n— l)£/ci. r indicates the duration 
in which the longitudinal sound wave propagates from 
one end of the chain to the other. Let us call the time 
t c during which the collision takes place "contact time". 
Here we note that contact time is almost equal to 2r 




Dimensionless time 

FIG. 2: Trajectories of the particles in the n — 20 colliding 
chain. (Reduced units where c; = r = 1.) 



for large n. COR and contact time arc independent of 
the collision velocity of the chain because time intervals 
At q = t q+ i — t q are independent of initial velocity for 
every q and the orbits retain similar forms even if initial 
velocity of the chain varies. Therefore, Ej/E, the ratio 
of the j-th fundamental mode energy to total energy, 
does not depend on initial velocity. The COR, which is a 
ratio of energy, also does not depend on initial velocity. 
Hereafter, we use the reduced unit c/ = r = 1, setting I, 
m and k as unity. 

In model A, the coefficient of restitution 77 depends 
only on n. An approximate expression of rj was derived 
by Basile and Dumont |l~5|] . To simplify the problem, we 
assume that At q and the velocity before the g-th colli- 
sion v q are constant for every q. Numerical simulations 
show that these are good approximations. Under this 
assumption, the vibration energy of the j-th mode after 
the collision process is given as 



E; 



Am 



2 ? 
V cos 



vrA sin 2 {^/(n)At/2} 
2n) sin 2 {wjAi/2} : 



(7) 



where uij is the frequency of the j-th mode, which is 
ujj = sin(7rj/n). At, fjn) and v are determined from 
numerical simulations Jig ] . However we can also estimate 
these values by solving the collision of the chain with 
n = 2: 



At = 



fin) 



V2u> 

2^2 
r 

7T 
7T 

2\/2 



2.22/w (2.31/w), (8) 



0.90n (0.867n), (9) 
«b = l.ll«D (1.15wq), (10) 



where the values in brackets were obtained from a nu- 
merical simulation of an n — 500 chain. We can estimate 
an approximate value of 77 by using Eqs. (||), (^|) and 
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FIG. 3: A log-log plot of dissipation rate versus number of 
particles. The solid line has slope —2/3. (Reduced units 
where c; = r = 1.) 



FIG. 4: Coefficient of restitution for the collisions of log-, 
Toda- and Lennard-Jones- type chains as a function of collision 
velocity. The chains consist of 100 particles. 



To obtain the asymptotic behavior of 1 — r/ for large n, 
we expand the dispersion relation. 



7TJ 7T 2 j 2 

m — 1 



(11) 



Substituting only the leading order of Eq.(^Tj) into 
Eq.(^), we have rj = 1. Taking the second order into 
account, we obtain 



1 - rj ~ 0.652 n~3, 



(12) 



in the limit n — > oo. This relationship agrees very well 
with the numerical result. In Fig. ||, we plot the dissipa- 
tion rate 1 — 77 versus n. 



III. UNHARMONIC CHAINS 

In the case of harmonic chains, COR and contact time 
were independent of initial velocity vq- Here, we intro- 
duce nonlinearity to the springs in the chain (model B). 

In order to maintain a universal viewpoint, we first 
consider a velocity scale characterized by the nonlinear- 
ity of the spring. Let U (x) be the potential energy of the 
spring, which is chosen to be U(l) — assuming the nat- 
ural length of the spring to be unity. We can define the 
amplitude x* at which the harmonic term and the sum of 
the remaining terms are equal in the Taylor expansion, 
i.e., x* is given by the solution of 



1 d 2 U{x) 
2 



dx 2 



x=l 



(x* - I)' 



The velocity which corresponds to x* is 



v* = y / 2U(x*). 

Hereafter, we discuss the velocity dependence of the col- 
lisions on the scale of v*. 



Let us choose the following three types of potential and 
compare the results of simulation, 
(a) Lennard-Jones potential 




(b) Toda potential 

U toda (x) = ft e -'(*-i) +0(3.-1) 


(c) Log-type potential 

Ulog{x) = X-\0g(x) 

For the Toda potential, we set ab = 1 and b — 10. Each 
potential has one minimum at x = 1 and the function 
forms are similar. However, increasing behaviors of the 
repulsive forces derived from the three potentials are, in 
a very short distance, different from one another. In Fig. 
||, we plot the COR versus initial velocity for each case, 
as determined by numerical simulations. In the limit of 
small velocity, it is clear that COR approaches the value 
obtained for the harmonic chain. The COR decreases 
with increasing initial velocity and has a minimum near 
v = 1. We note that the COR lies on almost the same 
curve for all three types of potential. 

Figure |5| shows n-dependence of r\ for different values of 
initial velocity. In very low-velocity collision, as expected 
from the result for model A, dissipation rate 1 — 77 will ap- 
proach to zero in the thermodynamic limit. Near vq = 1, 
the COR does not approach unity but remains a con- 
stant less than unity even in the thermodynamic limit. 
We can hence conclude that nonlinearity of the potential 
U (x) causes the dissipation, which does not appear in 
model A. 

Using a technique based on the perturbative theory, we 
consider the collision of model B for small initial incident 
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FIG. 5: The relationship between dissipation rate and number 
of particles in the collision of Lennard- Jones-type chain. Each 
plot corresponds to a different initial velocity: vg = 0.780 (a), 
0.390 (0), 0.160 (♦), 0.078 (□) and 0.008 (•). 



velocity » « 1. During collision processes, the particle 
j = 1 transmits vibration force to its neighbor parti- 
cle. Let us regard this force as external force F(t) which 
acts on the chain. The characteristic frequency of this 
force is io ex t — 2-ir/At ~ 2\f2tn. This frequency is higher 
than any frequency of fundamental modes of the chain. 
Hence, no fundamental mode is excited by the force. In 
this situation, the amplitude of each particle's vibration 
is rapidly damped progressively into the chain, i.e., the 
particle j = 1 has the largest amplitude during collision. 
It is expected that At will shorten with increasing initial 
velocity. As a first approximation, we take into account 
only the change of At against initial velocity vq in Eq. 
(R). To estimate the velocity dependence of At, let us 
consider the collision of a Lennard- Jones-type chain with 
n = 2. In this case, particle j = 1 collides with the wall 
two times. The Hamiltonian is 

n = ^(xj+x 2 2 ) + U L j{x 2 -x 1 -l) + V w . (13) 

Let x = X2 — x\ — i and Xq — {x% + x{)/2. For initial 
conditions, we adopt Eq. (g). Immediately after the first 
collision of particle x%, x and x g and their time deriva- 
tives become 





-2v , 



1/2 




(14) 



Taking into account the second order of Ulj and solving 
the equation of motion with the initial conditions of Eq. 
( |l4| ) in the first-order perturbation theory, one can ob- 
tain the interval At between first and second collisions 
of particle x\ as 



At: 



n \ v*j 



V2 



(15) 



FIG. 6: Increase of dissipation rates. We plot the dissipation 
rate of model B minus that of model A for r = 3, 4 and 5. 
Fitted lines have slopes 1.34,2.09 and 3.14, respectively. 



with a — 0.808. Substituting this into Eq. 



we have 



E; 



COS 



2 r l 



Wj/(n)At(l - av)] 



sin 



2 r l 



L2 w j At(l — av)] 



(16) 



In the limit n 
form 



oo, u>j can be replaced with the linear 

(17) 



When v/ci < 1, Eq. @ can be approximated as 



32 
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Substituting this into Eq. 
rate 



El 

E 



64 



64 

Cv, 



we have the dissipation 



dx—^ sin 2 (avnx) 



dx 



sin 2 (ax) 



(18) 



where C = -p-a ~ 2.61. This implies that the dissipa- 
tion rate 1 — r] increases with the power law v p where 
p = 1. This result agrees with numerical simulation for 
low-velocity collision. However the constant C is not in 
accord with the above result (our numerical simulation 
gives C — 0.66). 

The above result directly depends on the exponent of 
the lowest un-harmonic term of interparticle potential 
energy. In the case of Lennard- Jones potential, the ex- 
ponent r — 3. The exponent is the same for other types 
of chains. Let us discuss more general cases. Suppose 
that the interparticle potential energy can be expanded 
around its equilibrium position as 

U{x) = \x 2 + cx r + higher order, (19) 
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Dimensionless velocity 

FIG. 7: Plot of ratio of f(n) to n and COR versus velocity 
for Lennard- Jones- type chain. Both increase at v « 1. 



type. Collision between the two spheres is reduced to a 
simple exchange of their momenta. We can exactly solve 
the entire dynamics of the collision between the chain 
and a rigid wall. The sphere j = 1 collides with the 
wall n times and no internal vibration remains after the 
collision, i.e., the COR is exactly unity in this case. 

In the limit of high- velocity collision of model B, par- 
ticles interact like rigid spheres because Ulj, U to da and 
Ui g all behave like hard core potentials within a very 
short distance. In Fig. 0, we show the COR of a Lennard- 
Jones-type chain and f(n)/n versus collision velocity on 
the same plot. The rate f(n)/n approaches unity for 
v > 1. This indicates that particles interact like rigid 
spheres in high-velocity collision. Consequently COR in- 
creases in the high-velocity regime. This is a feature 
which only the one-dimensional model exhibits and is un- 
realistic. In real systems, plastic deformation is crucial 
in such high- velocity collisions. 



where the constant c is a positive (negative) number when 
r is even (odd). In this case, the contact time shortens 



as v 



The dissipation rate hence increases as 



r-2 



(20) 



Taking potential energy Eq.(|l9|) for the interaction of 
particles in model B, we plot numerical results of COR 
in Fig. ^ for v <C 1. The results agree with Eq. (20). 

When the particles in the chain can be regarded as 
rigid spheres, interparticle interaction is of the 5 function 



IV. SUMMARY 

We have presented a simple one-dimensional micro- 
scopic model of colliding bodies to understand the energy 
dissipation process. Lennard-Jones, Toda and Log-type 
potentials are chosen as interactions between particles. 
We found that the COR depends on the initial velocity 
and is minimum at v/v* ~ 1. These behaviors are inde- 
pendent of the potential form. In low- velocity collisions, 
the relationship between the energy dissipation rate and 
collision velocity is derived using perturbation methods. 
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